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We study the proximity effect between an s-wave superconductor and the surface states of a strong 
topological insulator. The resulting two dimensional state resembles a spinless p x + ip y supercon- 
ductor, but does not break time reversal symmetry. This state supports Majorana bound states 
at vortices. We show that linear junctions between superconductors mediated by the topological 
insulator form a non chiral 1 dimensional wire for Majorana fermions, and that circuits formed from 
these junctions provide a method for creating, manipulating and fusing Majorana bound states. 
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Excitations with non-Abelian statistics [l[ are the ba- 
sis for the intriguing proposal of topological quantum 
computation^. The simplest non-Abelian excitation is 
the zero energy Majorana bound state (MBS) associated 
with a vortex in a spinless p x +ip y superconductor [1,0, S 
6]. The presence of 2N vortices leads to a 2 N fold ground 
state degeneracy. Braiding processes, in which the vor- 
tices are adiabatically rearranged, perform non trivial 
operations in that degenerate space. Though MBSs do 
not have the structure necessary to construct a universal 
quantum computer[3], the quantum information encoded 
in their degenerate states is topologically protected from 
local sources of decoherencejij]. 

MBSs have been proposed to exist as quasiparticle ex- 
citations of the v = 5/2 quantum Hall effect [j], in 
the cores of h/4e vortices in the v- wave superconduc- 
tor Sr 2 Ru049] and in cold atoms 13, 11|. In this paper 
we show that the proximity effect between an ordinary 
s-wave superconductor and the surface of a strong topo- 
logical insulator (TI)[l2, 13, 14, 15| leads to a state which 
hosts MBSs at vortices. We then show that a linear 
superconductor - TI - superconductor (STIS) junction 
forms a non chiral ID wire for Majorana fermions. Such 
junctions can be combined into circuits, which allow for 
the creation, manipulation and fusion of MBSs. 

A strong TI is a material with an insulating time rever- 
sal invariant bandstructure for which strong spin orbit in- 
teractions lead to an inversion of the band gap at an odd 
number of time reversed pairs of points in the Brillouin 
zone. Candidate materials include the semiconducting 
alloy Bii-^Sbo;, as well as HgTe and a-Sn under uniaxial 
strain[l5l]. Strong TIs are distinguished from ordinary 
insulators by the presence of surface states, whose Fermi 
arc encloses an odd number of Dirac points and is asso- 
ciated with a Berry's phase of n. In the simplest case, 
there is a single non degenerate Fermi arc described by 
the time reversal invariant Hamiltonian 



Hq = ijy(—ivd ■ V — p)ip. 



(1) 

Here ijj = (tp^,^p{) T are electron field operators, a = 
(a x ,cr y ) are Pauli spin matrices and p is the chemical 



potential. Hq can only exist on a surface because it vio- 
lates the fermion doubling theorem[l6j]. The topological 
metal is essentially half of an ordinary 2D electron gas. 

Suppose that an s-wave superconductor is deposited 
on the surface. Due to the proximity effect, Cooper pairs 
can tunnel into the surface states. This can be described 
by adding V = A^jvj + h.c. to H Q , where A = A e^ 
depends on the phase 4> of the superconductor and the 
nature of the interface [171]. The states of the surface can 
then be described by H = ^H'$>/2, where in the Nambu 
notation = ((V>|, ipi), (V>|, — V'{)) T ar >d 

H = -ivT z o ■ V — \it z + A (t x cos(/> + t v sin</>). (2) 

r are Pauli matrices that mix the ip and tjy blocks of 
Time reversal invariance follows from [0,7i] = 0, 
where O = ia y K and K is complex conjugation. Par- 
ticle hole symmetry is expressed by S ~ a v r y K, which 
satisfies {E,H} = 0. When A is spatially homogeneous, 
the excitation spectrum is = ±\/ (±v\k\ — /i) 2 + A 2 ,. 
For p 3> Ao, the low energy spectrum resembles that 
of a spinless p x + ip y superconductor. This analogy can 
be made precise by defining = ("0Tk + e lSk, 0|k) / \/2 for 
k = fco(cos 6k, sin 9k) and vko ~ fi. The projected Hamil- 
tonian is then ^kM^I — A*) c jk c k + (Ae l6 ' k cj c cL k + /i.c.)/2. 
Though this is formally equivalant to a spinless p x + ip y 
superconductor there is an important difference: TL re- 
spects time reversal symmetry, while the p x + ip y super- 
conductor does not. 

It is well known that a h/2e vortex in a p x + ip y su- 
perconductor leads to a MBS [!]• This suggests that for 
/i ^> Ao a similar bound state should exist for ([2]). The 
bound states at a vortex are determined by solving the 
Bogoliubov de Gennes (BdG) equation H£, = E£ in polar 
coordinates with A(r, 0) — Ao(r)e ±l6 '. A zero energy so- 
lution exists for any p. The algebra is simplest for p = 0, 
where the zero mode has the form 



± - f dr'A (r')/v 



(3) 



with x+ = ((0, i), (1, 0)) T and x ~ = ((1,0), (0, ~i)) T . 

Another feature of p x +ip„ superconductors is the pres- 
ence of chiral edge states [3, 0, El] • With time rever- 
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FIG. 1: (a) A STIS line junction, (b) Spectrum of a line junc- 
tion for W = = as a function of momentum for various 
<f>. The solid line shows the Andreev bound states for <f> = it. 
The dashed lines are for <f> = 37r/4, it/2 and n/4. The bound 
states for 4> — merge with the continuum, indicated by the 
shaded region, (c) A tri-junction between three superconduc- 
tors, (d) Phase diagram for the tri-junction. In the shaded 
regions there is a ± MBS at the junction. 



sal symmetry, chiral edge states can not occur in our 
system. The surface - which itself is the boundary of 
a three dimensional crystal - can not have a boundary. 
By breaking time reversal symmetry, however, a Zccman 
field can introduce a mass term Ma z into I|ll2p which can 
open an insulating gap in the surface state spectrum. By 
solving J2]) we find that the interface between this in- 
sulating state and the superconducting state has chiral 
Majorana edge states. This could possibly be realized 
by depositing superconducting and insulating magnetic 
materials on the surface to form a superconductor-TI- 
magnet (STIM) junction. It is interesting to note that 
for spinless electrons the p x + ip y superconductor vio- 
lates time reversal, while the vacuum does not. For our 
surface states it is the insulator which violates time re- 
versal. A related effect could also occur at the edge of 
a two dimensional TI 2^, 21, 22|. which is described by 
(|1I2|) restricted to one spatial dimension. At the bound- 
ary between a region with superconducting gap At x and 
a region with insulating gap Ma z we find a MBS, anal- 
ogous to the end states discussed in Refs. |23|, [2J. In the 
following we will focus on STIS junctions, which can lead 
to non chiral one dimensional Majorana fermions, as well 
as MBSs. 

Consider a line junction of width W and length L — > oo 
between two superconductors with phases and <p in 
contact with TI surface states. We analyze the An- 
dreev bound states in the surface state channel between 
the superconductors by solving the BdG equation with 
A(x,y) = A e^ for y > W/2, A for y < -W/2 and 
otherwise. The calculation is similar to Titov, Ossipov 
and Beenakker's[25| analysis of graphene SNS junctions, 



except for the important difference that graphene has 
four independent Dirac points, while we only have one. 
For W <C v/Aq there are two branches of bound states, 
which disperse with the momentum q in the x direction. 
For W — ji — we find 



E±(q) = ± 



A 2 cos 2 (0/2)] 



1/2 



(4) 



For <fi = 7T the spectrum is gapless. It is useful to 
construct a low energy theory, for q ~ and 4> — n ~ e - 
Finite W and fi can then easily be included. We first solve 
the BdG equation for the two E = modes Ca=i,2(y) at 
q = and <p = it. It is useful to choose them to satisfy 
S£ a = £ a . Up to a normalization they may be written 



Cl±<2 = ((l,±i),(±Z,-l)) 5 



^±ifJ.y/v-J^ yl dyA (v)/v 



(5) 



We next evaluate {C,a\qv x T z \C,b) and (( a \eA 9(y 
W)r y \£b) to obtain the "k ■ p" Hamiltonian, 



n 



-ivT x d x + St v , 



(6) 



where v = v[cos fxW + (A //j,) sin/^F] A 2 ,/^ 2 + A 2 ,) and 
S = A cos(</>/2). The Pauli matrices act on ( a and 
are different from those in {2]). In this basis 8 = ir y K 
and S = K. 7i resembles the Su Schrieffer Heeger (SSH) 
model[26j]. However, unlike that model, the E±(q) states 
are not independent, and the corresponding Bogoliubov 
quasiparticle operators satisfy 7+ (q) = j-(—q)^. The 
system is thus half a regular ID Fermi gas, or a non 
chiral "Majorana quantum wire" . 

Below it will be useful to consider junctions that bend 
and close. When a line junction makes an angle 9 with 
the x axis the basis vectors (JSJ are modified according 



to ( a 



Jo-,,8/2 



Cq- H, however, is unchanged even when 



0(x) varies. On a circle, £ a changes sign when 9 advances 
by 2-7T. Therefore, eigenstates of TL must obey antiperiodic 
boundary conditions, ip(0) — —ip(2n). 

Next consider a tri-junction, where three superconduc- 
tors separated by line junctions meet at a point, as in Fig. 
lc. When ^=1,2 is m the shaded region of Fig. Id, a 
MBS exists at the junction. Though the general BdG 
equation cannot be solved analytically, this phase dia- 
gram can be deduced by solving special limits. When 
4>k = there is no bound state. Another solvable limit 
is when three line junctions with W — are oriented at 
120°, and (f>k = ±fc(27r/3). This is a discrete analog of a 
± vortex with C3 symmetry, and is indicated by the cir- 
cles in Fig. Id. For /i = we find a MBS identical to © 
with the exponent replaced by —Aqti ■ r/v. Here n is a 
constant unit vector in each superconductor that bisects 
the angle between neighboring junctions. The MBS can 
not disappear when (f>k are changed continuously unless 
the energy gap closes. The phase boundaries indicated 
in Fig. Id therefore follow from the solution of the line 
junction, and occur when the phase difference between 
neighboring superconductors is tt. 
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FIG. 2: Energy levels in units of AE = v/L for a STIS line 
junction terminated by two tri-junctions as a function of 8 = 
Ao cos 0/2 for ~ 7r. In (a) two MBSs are created or fused 
when passes through 7r. In (b) a single MBS is transported 
from one end to the other. The insets depict the MBSs. 



It is instructive to consider the limit where two of the 
lines entering the tri-j unction are nearly gapless. For 
(f>i = 7T — ei and 4>2 — 7T— £2 Fig- Id predicts a MBS when 
e±€2 < 0. This can be understood with Eq. [SJ which 
describes the lower two line junctions, which have masses 
(5i. 2 = A ei, 2 /2. When ei£2 < <5 changes sign, leading 
to the well known midgap state of the SSH model[26l.l27||. 
which in the present context is a MBS. 

A line junction terminated by two tri-junctions allows 
MBSs to be created, manipulated and fused. When 
passes through it MBSs appear or disappear at both 
ends. To model this we assume the phases of the super- 
conductors on either side of the line junction are and 
7T— e, and that the superconductors at the left (right) ends 
have phases <Pl(r), which are not close to or ir. This 
allows us to model the ends using a hard wall bound- 
ary condition S — > ±oo, where the sign at each end is 
s l,r — sgnsin^L^. It is straightforward to solve $Q to 
determine the spectrum as a function of 8 — Aoe/2 for 
a line of length L using this boundary condition. There 
are two cases depending on the sign of slSr. 

For sl = Sfi = 1 either zero or a i pair of MBSs 
are expected. The spectrum, shown in Fig. 2a, may be 
written E±(S) = ±yJW 
tan k n L 
For L - 



■ v 2 k 2 , where k n 



are solutions to 

; ~vk n /8. Midgap states are present for 8 < 0. 
oo a pair of zero energy states £1,2(2,2/) are 
localized at each end with wavefunctions 



6 = Cie 



6 = C2e" 



\S\(L-x)/v 



(7) 



where ( a (y) are given in (J5]). For finite —8L/v ^ 1 the 
eigenstates are ip± = £1 ± with energies Eq(8) ~ 
±2|(5|e~l d l L /' u . These define Bogoliubov quasiparticle op- 
erators, r± = ((p±)' i If. Since <p± = S<^ T it follows that 
r + = ri = (71 — «72)/2 where j a — (£ a )^ are Majorana 
operators. The pair 71,2 thus define a two state Hilbcrt 
space indexed by n.12 = r^_r + . The splitting between ip± 
then characterizes the interaction between the MBSs, 

H = E+(5)(Ylv + - 1/2) = ^ +(«5)7 27l /2. (8) 
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FIG. 3: Simple circuits made from STIS junctions. When <f> 
is advanced from to 2-7T (a) produces an entangled state and 
(b) interchanges MBSs 1 and 2. 



The Sl = sr = —1 case is similar. Eq. |8] applies to both 
cases, provided 72 is associated with the + vortex. 

This provides a method for both creating and fusing 
pairs of MBSs. Suppose we begin in the ground state 
at 8 > with no MBSs present. Upon adiabatically 
decreasing 8 through MBSs appear in the state |0i 2 ). 
Next suppose that initially 8 < 0, and a pair of MBSs 
are present in the state |ni2). When 8 is adiabatically in- 
creased through the system will remain in \ri12), which 
will either evolve to the ground state or to a state with 
one extra fermion. The difference between the two states 
can be probed by measuring the current flowing across 
the linear junction, which depends on whether the An- 
dreev bound state tp+ is occupied. The measured current 
will be / = J± A//2, where the current carried by (p+ 
is AI = {e/h)dE+ /d(j> - eA /27i for SL/v > 1. For 
A - .1 meV AI - 10 nA. 

Finally, consider the case sl = —sr = 1, in which 
a — MBS is at one end or the other, as in Fig. 2b. 
There are plane wave solutions with energy E^ — 
±y/S 2 + (nirv/L) 2 for n 



= 1,2,. 

Eq = state with wavefunction 

£(x,y) oc (i(y)e 



energy 
along with a single 



5x I v 



(9) 



Depending on the sign of 8, £ is exponentially localized 
at one end or the other. When 8 changes sign, the MBS 
smoothly switches sides. This provides a method for 
transporting a MBS from one node to another. 

We now discuss simple circuits built from STIS junc- 
tions. First, consider Fig. 3a and a process in which the 
phase of the central island is adiabatically advanced from 
to 2tt. For = there are no MBSs. At (j) = 2vr/3 
two pairs of MBSs are created at the top and bottom 
line junctions. At <f> = 4tt/3 the MBSs are fused at the 
left and right line junctions. If the system begins in the 
ground state = 0, then when — > 27r we findf^ 

IO12O34) -> (IO14O32) + |luls2»/V2. (10) 

Thus, after the cycle, the left and right segments are in 
an entangled state. The currents measured across the left 
and right junctions will be 7± AJ/2 with 50% probability 
and will be perfectly correlated. 

Eq. [TO] can be understood in two ways. First, the cy- 
cle effectively creates two pairs of MBS's, interchanges 



4 



a pair (say 2 and 4) and brings the pairs back together. 
As shown by Ivanov 0] , this corresponds to the operator 
-Pm = (l + 7274)/v / 2, which leads directly to ([TO]), Alter- 
natively, this result can be derived from (|6I8[) . From ©, 
the Hamiltonian for (f> < 2-7r/3 is -f/i oc 1(7172 + 7374)- 
For <j> > 47r/3 it becomes #2 oc 1(7174 — 7372)- Here the 
minus sign arises because, as explained after Eq. [6] the 
closed ID circuit must have antiperiodic boundary con- 
ditions. Thus, one of the line junctions must have a cut 
where the wavefunction changes sign. We chose the cut 
to be on the junction between 2 and 3. It is then straight- 
forward to express the groundstate of Hi in terms of the 
eigenstates of H 2 , which leads directly to (flTJ)) . 

Fig. 3b gives a geometry for interchanging MBSs with- 
out fusing them. For 4> — MBSs are located as shown. 
When 4> advances by 2tt the MBSs hop counterclockwise 
three times and are interchanged. Ivanov's rulesQ pre- 
dict 72 — ► 71, 71 — * —72- Again the minus sign can be 
understood in terms of the cut due to antiperiodic bound- 
ary conditions. One can imagine larger arrays, where this 
process performs elementary braiding operations. 

The experimental implementation of this proposal will 
require progress on many fronts. The first is to find a 
strong TI with a robust gap. Bii_ x Sb^ and strained 
HgTe can have gaps of order 30 meV |l5[. The next 
is to interface with an appropriate superconductor. Ao 
depends on the quality of the interface, Schottky barriers 
and the mismatch in the Fermi wavelengths 17( . If these 
can be optimized, Ao can be comparable to the gap of 
the bulk superconductor [28j |. 

The simplest experimental geometry would be to con- 
sider a single line junction with W < hv/ Ao. For Ao ~ 
.1 meV and hv ~ 1 eVA this could be achieved with 
W < 1/im. This should be similar to a graphene SNS 
junction (25)]. A signature of the Majorana character of 
the junction could be probed by measuring the ther- 
mal conductance along the channel for fcgT < Ao. For 
<f> = 7r the central charge c = 1/2 of the gapless Majo- 
rana modes leads to a quantized Landauer thermal con- 
ductance k — c(ir 2 /3)(k%/h)T. By constructing a pair 
of tri-junctions as in Fig. 2 the presence of MBSs can be 
controlled. It would then be interesting to perform tests 
of the non locality of MBSs envisioned in Refs. 24, 29|. 

Manipulating and fusing MBSs places more stringent 
requirements on the energy gaps. The junctions should 
be sufficiently short that AE = v/L > ksT, but suffi- 
ciently long that the MBSs are well localized. The good 
news is that AE varies as a power of L, while the MBS 
overlap is exponential, so at low temperature both crite- 
ria can be achieved. 

If the process of varying 6 to manipulate the MBSs is 
non adiabatic or AE < ksT then additional quasiparti- 
cles could be excited. If those quasiparticles escape and 
interact with other MBSs, then the state of the MBSs 
will be disturbed. However, if 8 <C Aq the excited quasi- 



particles will be confined to the segment in which they 
were created. If 6 is turned up so that ksT <C 5 <C A , 
and the system relaxes back to its ground state, then the 
state of the MBSs will remain intact. Thus, if there is 
sufficient dynamic range between fcgT and Ao, the sys- 
tem can tolerate these excitations. 
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